Molecular-dynamics of water transport through membranes - water from solvent to solute

An application of Molecular Dynamics computer simulation (MD) to the process of transport of water through a lipid bilayer membrane is described. The permeation process is far too slow to be modeled by straightforward MD. In stead the inverse of the permeability coefficient is expressed as an integral over a local permeation resistance, which itself is inversely proportional both to the local density and the local diffusion constant. These quantities axe recovered from MD simulations. The local density relates to a free energy profile, which is constructed by a combination of density determination, of the mean force on constrained molecules, and particle insertion. Thus a slow process can be accurately predicted from relatively short MD simulations.</p

Per la corretta attribuzione del "Romanzo delle donne contemporanee in Italia" (1863)

The use of free energy simulation techniques in the study of protein stability is critically evaluated. Results from two simulations of the thermostability mutation Asn218 to Ser218 in Subtilisin are presented. It is shown that components of the free energy change can be highly sensitive to the computational details of the simulation leading to the conclusion that free energy calculations cannot currently be used to reliably predict protein stability. The different factors that undermine the reliability are discussed

Multi-Scale Simulations Provide Supporting Evidence for the Hypothesis of Intramolecular Protein Translocation in GroEL/GroES Complexes

The biological function of chaperone complexes is to assist the folding of non-native proteins. The widely studied GroEL chaperonin is a double-barreled complex that can trap non-native proteins in one of its two barrels. The ATP-driven binding of a GroES cap then results in a major structural change of the chamber where the substrate is trapped and initiates a refolding attempt. The two barrels operate anti-synchronously. The central region between the two barrels contains a high concentration of disordered protein chains, the role of which was thus far unclear. In this work we report a combination of atomistic and coarse-grained simulations that probe the structure and dynamics of the equatorial region of the GroEL/GroES chaperonin complex. Surprisingly, our simulations show that the equatorial region provides a translocation channel that will block the passage of folded proteins but allows the passage of secondary units with the diameter of an alpha-helix. We compute the free-energy barrier that has to be overcome during translocation and find that it can easily be crossed under the influence of thermal fluctuations. Hence, strongly non-native proteins can be squeezed like toothpaste from one barrel to the next where they will refold. Proteins that are already fairly close to the native state will not translocate but can refold in the chamber where they were trapped. Several experimental results are compatible with this scenario, and in the case of the experiments of Martin and Hartl, intra chaperonin translocation could explain why under physiological crowding conditions the chaperonin does not release the substrate protein

Long range molecular dynamics study of regulation of eukaryotic glucosamine-6-phosphate synthase activity by UDP-GlcNAc

Glucosamine-6-phosphate (GlcN-6-P) synthase catalyses the first and practically irreversible step in hexosamine metabolism. The final product of this pathway, uridine 5’ diphospho N-acetyl-D-glucosamine (UDP-GlcNAc), is an essential substrate for assembly of bacterial and fungal cell walls. Moreover, the enzyme is involved in phenomenon of hexosamine induced insulin resistance in type II diabetes, which makes it a potential target for antifungal, antibacterial and antidiabetic therapy. The crystal structure of the isomerase domain of GlcN-6-P synthase from human pathogenic fungus Candida albicans, in complex with UDP-GlcNAc has been solved recently but it has not revealed the molecular mechanism of inhibition taking place under UDP-GlcNAc influence, the unique feature of the eukaryotic enzyme. UDP-GlcNAc is a physiological inhibitor of GlcN-6-P synthase, binding about 1 nm away from the active site of the enzyme. In the present work, comparative molecular dynamics simulations of the free and UDP-GlcNAc-bounded structures of GlcN-6-P synthase have been performed. The aim was to complete static X-ray structural data and detect possible changes in the dynamics of the two structures. Results of the simulation studies demonstrated higher mobility of the free structure when compared to the liganded one. Several amino acid residues were identified, flexibility of which is strongly affected upon UDP-GlcNAc binding. Importantly, the most fixed residues are those related to the inhibitor binding process and to the catalytic reaction. The obtained results constitute an important step toward understanding of mechanism of GlcN-6-P synthase inhibition by UDP-GlcNAc molecule

Molecular dynamics simulations reveal that AEDANS is an inert fluorescent probe for the study of membrane proteins

Computer simulations were carried out of a number of AEDANS-labeled single cysteine mutants of a small reference membrane protein, M13 major coat protein, covering 60% of its primary sequence. M13 major coat protein is a single membrane-spanning, α-helical membrane protein with a relatively large water-exposed region in the N-terminus. In 10-ns molecular dynamics simulations, we analyze the behavior of the AEDANS label and the native tryptophan, which were used as acceptor and donor in previous FRET experiments. The results indicate that AEDANS is a relatively inert environmental probe that can move unhindered through the lipid membrane when attached to a membrane protein

Detailed Regulatory Mechanism of the Interaction between ZO-1 PDZ2 and Connexin43 Revealed by MD Simulations

The gap junction protein connexin43 (Cx43) binds to the second PDZ domain of Zonula occludens-1 (ZO-1) through its C-terminal tail, mediating the regulation of gap junction plaque size and dynamics. Biochemical study demonstrated that the very C-terminal 12 residues of Cx43 are necessary and sufficient for ZO-1 PDZ2 binding and phosphorylation at residues Ser (-9) and Ser (-10) of the peptide can disrupt the association. However, only a crystal structure of ZO-1 PDZ2 in complex with a shorter 9 aa peptide of connexin43 was solved experimentally. Here, the interactions between ZO-1 PDZ2 and the short, long and phosphorylated Cx43 peptides were studied using molecular dynamics (MD) simulations and free energy calculation. The short peptide bound to PDZ2 exhibits large structural variations, while the extension of three upstream residues stabilizes the peptide conformation and enhanced the interaction. Phosphorylation at Ser(-9) significantly weakens the binding and results in conformational flexibility of the peptide. Glu210 of ZO-1 PDZ2 was found to be a key regulatory point in Cx43 binding and phosphorylation induced dissociation

A Structural Model of the Pore-Forming Region of the Skeletal Muscle Ryanodine Receptor (RyR1)

Ryanodine receptors (RyRs) are ion channels that regulate muscle contraction by releasing calcium ions from intracellular stores into the cytoplasm. Mutations in skeletal muscle RyR (RyR1) give rise to congenital diseases such as central core disease. The absence of high-resolution structures of RyR1 has limited our understanding of channel function and disease mechanisms at the molecular level. Here, we report a structural model of the pore-forming region of RyR1. Molecular dynamics simulations show high ion binding to putative pore residues D4899, E4900, D4938, and D4945, which are experimentally known to be critical for channel conductance and selectivity. We also observe preferential localization of Ca2+ over K+ in the selectivity filter of RyR1. Simulations of RyR1-D4899Q mutant show a loss of preference to Ca2+ in the selectivity filter as seen experimentally. Electrophysiological experiments on a central core disease mutant, RyR1-G4898R, show constitutively open channels that conduct K+ but not Ca2+. Our simulations with G4898R likewise show a decrease in the preference of Ca2+ over K+ in the selectivity filter. Together, the computational and experimental results shed light on ion conductance and selectivity of RyR1 at an atomistic level

Configurational Entropy and Diffusivity of Supercooled Water

We calculate the configurational entropy S_conf for the SPC/E model of water for state points covering a large region of the (T,rho) plane. We find that (i) the (T,rho) dependence of S_conf correlates with the diffusion constant and (ii) that the line of maxima in S_conf tracks the line of density maxima. Our simulation data indicate that the dynamics are strongly influenced by S_conf even above the mode-coupling temperature T_MCT(rho).Comment: Significant update of reference

Visualization of Early Events in Acetic Acid Denaturation of HIV-1 Protease: A Molecular Dynamics Study

Protein denaturation plays a crucial role in cellular processes. In this study, denaturation of HIV-1 Protease (PR) was investigated by all-atom MD simulations in explicit solvent. The PR dimer and monomer were simulated separately in 9 M acetic acid (9 M AcOH) solution and water to study the denaturation process of PR in acetic acid environment. Direct visualization of the denaturation dynamics that is readily available from such simulations has been presented. Our simulations in 9 M AcOH reveal that the PR denaturation begins by separation of dimer into intact monomers and it is only after this separation that the monomer units start denaturing. The denaturation of the monomers is flagged off by the loss of crucial interactions between the α-helix at C-terminal and surrounding β-strands. This causes the structure to transit from the equilibrium dynamics to random non-equilibrating dynamics. Residence time calculations indicate that denaturation occurs via direct interaction of the acetic acid molecules with certain regions of the protein in 9 M AcOH. All these observations have helped to decipher a picture of the early events in acetic acid denaturation of PR and have illustrated that the α-helix and the β-sheet at the C-terminus of a native and functional PR dimer should maintain both the stability and the function of the enzyme and thus present newer targets for blocking PR function

Generic mechanism for generating a liquid-liquid phase transition

Recent experimental results indicate that phosphorus, a single-component system, can have two liquid phases: a high-density liquid (HDL) and a low-density liquid (LDL) phase. A first-order transition between two liquids of different densities is consistent with experimental data for a variety of materials, including single-component systems such as water, silica and carbon. Molecular dynamics simulations of very specific models for supercooled water, liquid carbon and supercooled silica, predict a LDL-HDL critical point, but a coherent and general interpretation of the LDL-HDL transition is lacking. Here we show that the presence of a LDL and a HDL can be directly related to an interaction potential with an attractive part and two characteristic short-range repulsive distances. This kind of interaction is common to other single-component materials in the liquid state (in particular liquid metals), and such potentials are often used to decribe systems that exhibit a density anomaly. However, our results show that the LDL and HDL phases can occur in systems with no density anomaly. Our results therefore present an experimental challenge to uncover a liquid-liquid transition in systems like liquid metals, regardless of the presence of the density anomaly.Comment: 5 pages, 3 ps Fig